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CRITICAL STRESS OF THIN-WALLED CYLINDERS IN AXIAL COAIPRESSION 


By S. B. Batdorf, Ml"rrt ScHiLDCRonx, and Makxiel Steik 


SUMMARY 

Empirical design curves are presented jor the critical stress qf 
thin-u'alled cylinders loaded in axial compression. These 
curves are plotted in terms of the nondirnensional parameters qf 
small-deflection theory and are compared vnih theoretical curves 
derived for the buckling of cylinders unth simply supported and 
clamped edges. An empirical equation is given for the buckling 
of cylinders having a length-radius ratio greater than about 0.75. 

The test data obtained from various sources follow the general 
trend of the theoretical curve for cylinders with clamped edges, 
agreeing closely unth the theory in the case of short cylinders, but 
falling considerably below the theoretical results for long 
cylinders. The discrepancy in the case of long cylinders increases 
■with increasing values of the ratio of radius to wall thickness. 
Plotting curves for different values of this ratio reduces the 
scatter in the test data and a certain degree of correlation unth 
theory is achieved. Advantage is taken of this correlation to 
obtain estimated design curves for cylinders with simply sup- 
ported edges, for which little experimental information is 
available. 

REVIEW OF PREVIOUS WORK ON PROBLEM 

Solutions to the problem of the deterniination of the critical 
stress of thin-tvalled cylinders subjected to axial compression 
have been presented by a large number of authors. South- 
weh, Timoshenlio, Fliigge, and numerous other authors 
have obtained theoretical solutions by the use of the small- 
deSection theory. (See, for example, references 1 to 4.) The 
value given by the small-deflection theory for the buckling 
stress of a thin-waUed cyhnder of moderate length having 
simply supported edges is 


predicted theoretically; moreover, the observed buckle pat- 
tern is different from that predicted on the basis of theory. 
A number of attempts have been made to e.xplain these dis- 
crepancies theoretically. Fliigge (reference 3) considered the 
deviation of the actual edge supports from the support con- 
ditions assumed in the theoretical treatment. DonneU 
(reference 5) and also Fliigge considered the initial deviation 
from the perfect cylindrical shape. Neither of the two 
attempted explanations satisfactorily accounts for the dis- 
crepancy e-xisting between the theoretical and e.xperimenta] 
values of the buckling stresses of cylinders. 

Von Karmen and Tsien (reference 11) introduced a large- 
deflection theory to accoimt for the buckling behavior of 
long cylinders. They showed that a long cylinder can be in 
equilibrium in a buckled state at a stress that is much smaller 
than the critical stress of small-deflection theory and also 
succeeded in accounting for the buckle pattern observed in 
the early stages of buckling. Reference 11 suggested that 
when a cyhnder has an initial imperfection or is subjected to 
a shock, it might pass into one of these buckled states with- 
out ever having reached the critical load given by equa- 
tion (1). Based on the same approach, a theory for the 
budding stresses of perfect cylinders was proposed by Tsien 
(reference 12), which gave for loading by rigid screw-power 
testing machine 

Ff 

0.370 y (2) 

and for loading by ideal hydraulic testing machine or dead 
weight 

Et 

<r„= 0.238 y (.3) 


Et 




Et 

= 0.608— 


( 1 ) 


where 

Ox critical compressive stress 
E Young's modulus 
t wall thickness of cylinder 
r radius of cylinder 

p, Poisson's ratio (in the present paper p is taken to be 
0.316 whenever a value is assigned to it) 


The large-deflection theories fail in two respects to describe 
adequately the buckling behavior of actual cylinders. First, 
the theories are formulated only for long cylinders; equations 
(2) and (3) seriously underestimate the critical stress of very 
short cylinders. Second, even for long cylinders, attempts to 
determine e.xperimentahy the numerical coefficient C in the 
buckling formula 

(4) 


Experiments (references 6 to 10) have shown that the 
actual critical stress is much lower than that predicted by 
equation (1). Except in the case of short cyhnders, the 
e.xperiments usually give values only 15 to 50 percent of that 


have resulted in appreciable experimental scatter. The 
experimental scatter is due at least in part to the initial imper- 
fections of construction always present in real cylinders. 
(See fig. 1.) 


.54.3 
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'"Perfecf" cylinders 



Figure l.— Effect of initial defects upon tlie maximum compressive loads of struts and cylinders. (From reference 13. Ilg. 18.) 


In the absence of a complete and satisfactory theoretical 
solution for . the critical stress of cylinders, a number of 
authors have proposed empirical formulas derived from test 
data (references 6 to 8). One such formula, which takes into 
account the length of the cylinder, is due to Ballerstedt and 
Wagner (reference 8): . 


The first 


g-3.3(0 + O,c(I) 

parameter in this equation is appropriate 


for 


t . 


flat sheet and the second parameter — is included to take 


into account the effect of curvature. More recently 
Kanemitsu and Nojima (reference 9) compiled all available 
previous experimental results and conducted a number , of 
tests of their own. The formula of ^Yagner and Ballerstedt 
was modified in reference 9 to bring it into better agreement 
with experiment as follows: 



0.16 




( 6 ) 


Within its range of application ^0.1< — <1.5; 
500 g j ^3000^ equation (6) is in considerably better agree- 
ment with experiment than equation (5) but, because of 
the change in the exponents of the parameters, equation (6) 
does not have any rational basis and must be regarded as 


purely empirical. A complete divorce of theory and ex- 
periment, however, cannot be regarded as a satisfactory per- 
manent settlement of the problem, and the present report 
attempts to bring theory and experiment into reasonable 
accord. 

CONTRIBUTION OF PRESENT PAPER 


In the present paper the available test data for critical 
stresses "of cydinders are reexamiiu'd and theoretical results 
are used as a guide in fairing the curves, in extending the 
range of validity of the existing empirical results, and in 
achieving a more rational interpretation of tlu' test data, 
For this purpose the test data are plotted in terms of the 
parameters of cylinder theory^ and are compared with theo- 
retical results derived in the appendix on the basis of srnall- 
deflection theoiy. 

The cylinder-theory parameters used are 




<tJP 


and 

where 

/ El^ \ 

D flexural stiffness of plate per unit length 

L length of cylinder 
Z curvature parameter 

kx critical stress coefficient appearing in the equation 

kxir^D 
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The experimental data are used as the principal guide in 
determining the critical compressive stresses of long cylinders 
(large values of Z) and the theoretical results are used mainly 
to supplement the test data in determining the critical stress 
of very short cylinders (small values of Z). The experimental 
scatter is reduced by presenting different curves for cylinders 
with different values of the ratio of radius to waU thickness 
on the assumption that for long cylinders this ratio furnishes 
some indication of the initial imperfections of the cylinder. 
Although these curves were determined partly on the basis of 
theoretical considerations, they are for convenience referred 
to herein as empirical crnwes. 

RESULTS AND CONCLUSIONS 

The critical compressive stress for cylinders is given by 
the equation 


where the values of may be obtained from figure 2 for 
cylinders with either clamped or simply supported edges. 
The design curves for cylinders with clamped edges are 
established by the test results reported in references 5 to 9. 
(See fig. 3.) Each curve was faired through a series of 
test points w'hich were plotted for cylinders with nearly 
the same ratio of radius to wall thickness rjt. The estimated 
(dashed) parts of the design curves for simple support wern 
obtained by fairing between the known experimental curves 
for long cylinders (large values of Z), which according to 
theory should be the same whether the cylinders have 
simply supported or clamped edges, and the theoretical 
curves for very short cylinders (small values of Z ) . 



FiGf.'iiE 2.— Critical stress coefficients for thin-walled circular cylinders subjected to axial compression. 
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For long cylinders the buckling stress is considerably 
below the theoretical buckling stress, the amount of the 
discrepancy depending on the ratio of radius to wall thickness. 
For very short cylinders the values of the critical 
stresses approach those for flat plates (simply supported 
ends, ^ 1 = 1 ; clamped ends, icj:=4), for which the agreement 
between theoretical and experimental results is known to be 
good. The general trend of each empirical curve is similar 
to tliat of the theoretical curve, indicating the existence of a 
certain degree of correlation between theory and test data. 

At large values of Z, the curves for become straight 
lines given by the formula 

^"^= 1 . 15^2 ( 8 ) 

where C depends on the ratio of radius to wall thickness of 
the cylinders in the manner shown in figure 4. From 
equations (7) and (8) the following expression for the 


critical stress is obtained 

(7r= CE (9) 

Equations (8) and (9) may be used when the length of the 
cylinder is more than about of the radius. The empirical 
curves of reference 10 indicate that the critical stress is sub- 
stantially independent of length when the length is greatc*r 
than about % of the radius. (This result may be checked 

by noting that for Z>0.5 -j the experimental curves of 

figure 2 are substantially straight lines of unit slope.) 

In figure 5, the empirical formula of Kanemitsu and Nojima 
(equation (6) of the present paper, the best previously 
published formula for the buckling of cylinders) is plotted 
in terms, of the parameters and Z, The curves are cut 
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off at those values of Z corresponding to the lower limits of 
the range of dimensions within which the formula was in- 
tended to apply. In general, for the range covered, the 
curves are in reasonable agreement with the test data and 
with the curves of the present paper for cylinders with 
clamped edges. The practical importance of the present 
approach lies in the fact that the use of the theoretical 
parameters and the theoretical solutions as a guide in fairing 
the curves permits the removal of the lower limits on these 
curves and also permits estimated curves to.. be drayTi for 
the buckling stresses for simply supported cylinders, although 
experimental data are available only for cylinders with 
clamped edges. 

. . » f 

Langley Meiiorial Aehonautic.al Laboratory, 
Nation-AL Advisory Committee for Aeron.autics, 
Langley Field, Va., March 20, 1947. 






APPENDIX 

THEORY FOR CYLINDERS BUCKLING UNDER AXIAL COMPRESSION 


SYMBOLS 

m positive integer 

r radius of cylinder 

t wall thickness of cylinder 

w radial component of displacement, positive outward 
X axial coordinate of cylinder 

y circumferential coordinate of cylinder 

C coefficient appearing in a^=CE ^ 

/ \ 

D flexural stiffness of plate per unit length 

E Young’s modulus 

L length of cylinder 

Q operator on w defined in appendix 

Z curvature parameter T Vl— 

Um coefficient of deflection function 

ki critical-compressive-stress coefficient appearing in the 

k,T^ 


LH 

12Z^(m— 1)^ 


■{m—iyki 


formula cr*= 

Vm deflection function defined in the appendix 

«=^ 

X half wave length of buckles in circumferential direction 
n Poisson’s ratio 

(Tj critical axial compressive stress _ _ 
c>x^'by^ 

bx‘^^ “^dx^dy^' by* 
the inverse of V* defined by 

THEORETICAL SOLUTION 

The critical compressive stress at which buckling occurs 
in a cylindrical shell may be obtained by solving the equa- 
tion of equilibrium. 

Equation of equilibrium. — The equation of equilibrium for 
a slightly budded cylindrical shell under axial compression 
is (reference 14) 


I Et„_,b*w . ,d% 

-DVw-l-^jV + — 0 


(Al) 


where x is the coordinate in the axial direction and y is the 
coordinate in the circumferential direction. Dividing 
through equation (Al) by D gives 


12Z^ _,c>% , 


b^w 




(A2) 


where the dimensionless parameters Z and /lj- are defined by 

, (tJE 

The equation of equilibrium may be represented by 

Qw=0 (A3) 


where Q is defined by 


12ZA 




^ , 7^ bl 

+ 12 ^^2 


Method of solution. — Equation (A2) may be solved by use 
of the Galerkin method as outlined in reference 15. When 
this method is applied, the deflection to is expressed in series 
form as follows 

'y j O.TjiV'm 


w= 


7n*=I 


CA4) 


The set of functions Vm are chosen to satisfy the boundary 
conditions but need not satisfy the equation of equilibrium. 
The coefficients are determined by the equations 


nL 

VmQ.wdxdy={) 

Jo Jo 


(m=l,2, 3, . ..i) (A5) 


In the present paper the deflection functions wore chosen 
to satisfy the following conditions on w at the ends of the 
cylinder: For simply supported edges 


For clamped edges 


b^w 


bw 


Simply supported edges. — An expansion for w that is 
sinusoidal in the circumferential direction and perfectly 
general (subject to the boundary conditions for simple 
support) in the axial direction is 


to— 


■ try ^ . mwx 

=sm ^ ^am sm — 


(Afi) 
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where X is the half wave length of the buckles in the circum- 
ferential direction. (Equation (A6) is equivalent to equa- 
tion (A4) if 


.TV . TTItZ 
1 ;„=sm y sm^ 


(A7) 


W = SU1 y yOS COS 1 (A9) 

Each term of this series satisfies the conditions on w at the 
edges. The function E* is now defined as 


Substitution of e.xpressions (A6) and (AT) into equation (A5) 
and integration over the limits indicated give 

1 ^7-m^ 

(^=E 2, 3, . . .) (AS) 

where 


The minimum value of for a given Z is found by assum- 
ing a value for m and minimizing k^ with respect to This 
procedure is followed for various values of m until a minimum 
ki is reached. Figure 2 presents the theoretical critical 
stress coefficients for cylinders with simply supported edges 
subjected to axial compression. 

Clamped edges. — A procedure similar to that used for 
cylinders with simply supported edges may be followed for 
cylinders with clamped edges. The deflection function used 
is the foUowmg series 


=sm 




cos 


(m—l)TX 

L 


-cos 


(m+l)xy\ 


(rn- 


(AlO) 
= 1,2,3. . . .} 


After the same operations are carried out for clamped 
edges as those carried out for the case of simply supported 
edges, the foUowtng equations result ; 


For m=l, 


For m=2, 


ni(23A"T J-A) ~ OailA — -0 
Q2(ilAFdA) — 0 


For m=3, 4, 5, ... , 


'i 






(All) 


where 


— [ (m — 1 ) ^ + jS^H 


12Z^(m-l)* 

E[(m-l)^+d-? 


— (m — 1)- k^ 
(m—l, 2, 3, 


. -> 


ese equations have 

a solution if 

the following infinite determinant vanishes: 



m=l 


0 

-44 

0 

0 

0 

0 

0 

m=2 

0 

TA+dA 

0 

-.14 

0 

0 

0 

0 

m=3 

-fi/a 

0 

-14+44 

0 

-44 

0 

0 

0 

m=4 

0 

-44 

0 

44+44 

0 

-44 

0 

0 

m=5 

0 

0 

-44 

0 

44+44 

0 

-44 

0 

m=6 

0 

0 

0 

-.44 

0 

44+44 

0 

-.44 

rn = 7 

0 

0 

0 

0 

-44 

0 

44+44 

0 

m=8 

0 

0 

0 

0 

0 

-.44 

0 

-44 + -44q 


(A12) 


If the rows and columns are rearranged the infinite determinant can be factored into the product of two infinite subdeter- 
minants. The critical stresses may then be obtained from the following equation: 


m=l 

244+44 

-44 

0 

0 ... 

0 

0 

0 

0 ... 

m=3 

-x44 

44+44 

-44 

0 ... 

0 

0 

0 

0 ... 

m=5 

0 

-44 

44+-44 

-44 . . . 

0 

0 

0 

0 ... 

m=7 

0 

0 

-44 

44+44 . . . 

0 

0 

0 

0 ... 

m=2 

0 

0 

0 

0 ... 

-44+44 

-44 

0 

0 ... 

m=4 

0 

0 

0 

0 ... 

-44 

44+44 

-.44 

0 ... 

m— 6 

0 

0 

0 

0 ... 

0 

-44 

44+44 

-44 . . . 

m=8 

0 

0 

0 

0 ... 

0 

0 

-44 

44+-44q • • • 

• 

. 

. 

• 

. . . - 




. . . . 


= 0 (A13) 
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The infinite subdeterminant involving terms with odd 
subscripts corresponds to a symmetrical buckling pattern 
(a buckling pattern symmetrical about the plane perpendicu- 
lar to and bisecting the axis of the cylinder). The infinite 
subdeterminant involving terms with even subscripts cor- 
responds to an antisymmetrical buckling pattern. For 
brevity these subdeterminants will be referred to as the odd 
determinant and the even determinant, respectively. 

The first approximation is 


Odd determinant: 


25/i-[-jM3=0 

(A14) 

Even determinant: 


il/j-r A/4=0 

(A15) 

Tlie second approximation is 
Odd determinant: 



(Alb) 

Even determinant: 



(A17) 


These equations show that for a selected value of the curva- 
ture parameter Z the critical buckling stress of a cylinder 
depends upon the circumferential wave length. Since a 
structure bucldes at the lowest stress at which instability 
can occur, is minimized with respect to the wave length 
by substituting values of /3 into the equations until the mini- 
mum hx, can be obtained from a plot of /r* against /3. For a 
given Z the lower of the two values obtained from equations 
(Al4) and {A16) is the first approximation of the critical 
buckling stress and, similarly, the lower of the two. values of 
kx obtained from equations (Alb) and (A17) is the second 
approximation of the critical buckling stress. 

Figure 2 presents the theoretical critical stress coefficients 
for cylinders with clamped edges in axial compression as ob- 
tained from the second approximation, together with the 
exact curve for the case of simply supported edges. Although 
this solution is an upper-limit solution, the second approxi- 
mation for the critical stress coefficient of a cylinder with 
clamped edges must be very close to being exact" for inter- 
mediate and large values of Z because it is almost identical 
with the exact solution for a cylinder with simply supported 


edges, and the critical stress of a cylinder with clamped edges 
cannot be less than the critical stress for a cylinder with 
simply supported edges. For values of Z approaching 
zero, the accuracy of the second approximation is indicated 
by the fact that it coincides with the known exact solution 
(kx=^) for a long flat plate with clamped edges. 
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